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Abstract- 


I.  Introduction 

In  the  frarrework  of  the  construction  of  a  EEG-based  Brain 
Computer  I  nterface  (BCI )  it  was  suggested  that  EEC  patterns 
can  be  better  detected  with  EEC  data  transformed  with  the 
Surface  Laplacian  computation  (SL)  than  with  the 
unprocessed  raw  potentials  [1].  Flowe/er,  accurate  SL 
esti  rretes  requi  re  the  use  of  many  EEC  electrodes,  when  local 
estirration  methods  are  used  [2,3],  These  local  estirration 
methods  corrpute  the  SL  at  a  certai  n  electrode  position  on  the 
base  of  the  value  of  the  surrounding  nearest  electrodes.  This 
cause  arors  in  the  SL  estimation  at  the  electrodes  placed  at 
the  boundary  of  the  electrode  grid,  sirce  their  neighbors  are 
not  well  defined  [4,5].  The  requi rerrent  of  an  high  number  of 
electrodes  for  an  affordable  estirrate  of  the  SL  of  the  EEC 
distribution  would  pre/ent  the  application  of  SL  estimates  in 
a  BCI  de/ice  to  be  used  by  laypersons  in  real-life  conditions, 
due  to  practical  reasons,  namely  the  tirre  consuming 
procedure  of  the  scalp  electrode  positioning.  On  the  other 
hand,  there  ©cists  a  class  of  estimators  of  the  surface 
Laplacian  of  the  EEC  potential  distributions  that  is  based  on 
the  use  of  a  i  global!  computational  scherr©  such  as  the 
spherical  splines  [6,7],  in  which  the  surface  Laplacian  at  a 
certain  electrode  position  depends  from  the  values  at  all  the 
other  positions  of  the  recording  array.  From  these 
considerations,  in  this  paper  we  i  nvesti  gate  the  performaxes 
of  global  computational  methods  for  the  esti  mation  of  the  SL 
from  a  limited  number  of  electrodes,  based  on  the  spherical 
spline  approach.  The  working  hypothesis  at  the  base  of  the 


present  work  are  i)  the  use  of  global  interpolation  SL 
estimates  can  produce  reasonable  SL  waveforms  even  if  a 
reduced  number  of  electrodes  are  used  (low  resolution  SL); 
ii)  the  use  of  low  resolution  SL  waveforms  for  BCIs  allows 
percentage  of  classifications  of  mental  patterns  statistically 
similar  to  those  obtained  with  SL  waveforms  computed  with 
the  f  ul  I  recordi  ng  array.  I  n  order  to  i  nvesti  gate  the  f  i  rst  i  ssu© 
the  1  gold  standard!  of  the  SL  esti  mates  obtained  with  the  full 
recording  array  were  compared  to  the  spline- based  SL 
estimations  oby  ned  on  the  same  EEG  recordings  but  using 
only  9  electrodes,  uniformly  disposed  along  the  scalp  in  the 
position  of  the  international  1020  system  The  second  issue 
was  addressed  by  compari ng  the  recognition  rates  of  the  BCI 
system  obtained  using  low- resolution  SL  wavd'orms  with 
those  obtained  by  the  full  resolution  SL.  Flere,  we  report 
results  in  the  recognition  of  mental  patterns  with  two  linear 
classifiers  based  on  the  Signal  Space  Projection  (SSP) 
algorithm  [8]  and  Fisheris  linear  discriminant  functions 
[9,10].  The  interest  in  the  use  of  such  linear  classifiers  for 
BCI  i s  due  to  their  si rrple  training  and  decision  procedures. 
In  fact  these  procedures  do  not  involve  non-linear 
minimization  procedures  such  as  those  necessary  for  the 
neural  network  classifiers  already  used  in  the  BCI  field  [11- 
13].  This  of  course  was  at  the  expense  of  the  possibility  to 
separate  the  input  space  with  non  linear  discriminant 
functions.  Recognition  perforrrBnces  of  the  two  linear 
classifiers  on  unprocessed  and  SL-transforrred  EEG  data 
were  computed  from  a  group  of  five  healthy  people 
performing  two  rrotor- related  mertal  tasks,  narrely  imagined 
right  and  left  hand  movements. 


II.  Methodology 

A.  Data  Collection 

Five  healthy  subjects  (three  rreles  and  two  females) 
participated  voluntarily  in  experimants  where  th^  performed 
different  tasks,  i  ncl  udi  ng  the  i  magi  nati  on  of  the  movemert  of 
the  right  rriddle  finger  (Rl)  as  well  as  the  Id't  rriddle  finger 
(LI).  The  whole  scalp  was  covered  with  26  EEG  electrodes 
placed  onto  standard  locations  according  to  the  extension  of 
the  1020  international  system  Sampling  frequency  was 
400  FIz,  and  signal  was  bandpass  filtered  between  0.1  and 
100  FIz  before  digitization.  At  the  beginning  of  a  recording 
session,  subjects  remained  in  a  resting  stated  relax  with  ^es 
openedd  for  60s  The  EEG  activity  of  this  period  is  used  as  a 
b^ine  for  subsequent  analysis  of  the  mental  tasks  Then, 
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subjects  started  performing  the  task  immediately  after  the 
operator  instructed  them  to  do  so,  and  th^  maintained  that 
task  for  more  than  10s  Every  subject  executed  four  times 
each  task  during  the  recording  session,  with  a  resting  period 
of  10s  between  each.  After  removal  of  time  segments 
contaminated  with  EMG  in  the  arms  it  remains  about  40 
seconds  of  E  EG  si  gnal  s  for  each  task  for  e/ery  subj  ect. 

B.  SL  estimations. 

Surface  Laplacian  computations  were  performed  by  using  the 
spherical  spl  i  nes  of  order  two,  that  were  found  to  be  adequate 
to  describe  the  SL  distributions  of  EEG  data  [14],  For  each 
EEG  recording  two  computations  of  the  SL  were  perforrred. 
The  first  by  using  all  the  information  of  the  electrode 
montage  (full -resolution  SL),  while  the  second  SL 
computation  was  perforrred  by  using  only  the  data  of  the 
EEG  potentials  from  nine  electrodes  (low-resolution  SL). 
These  nine  electrodes  were  placed  on  the  scalp  surface 
according  to  the  position  F3,  Fz,  F4,  C3,  Cz,  C4,  P3,  Pz,  P4 
of  the  I  nternational  1020  System 

C.  Comparisons  of  the  SL  estimations 

The  aim  of  the  work  was  to  investigate  if  the  low- resolution 
SL  estirration  is  able  to  produce  relatively  accurate  SL 
estirretes  for  the  recognition  of  mental  patterns  Flence,  we 
used  the  correl  ati  on  cod'f  i  ci  ent  to  measure  the  f  i  t  between  the 
full  resolution  and  the  low-resolution  SL  waveforms  in  the 
nine  electrodes  selected.  These  comparisons  were  performed 
for  all  the  EEG  recordings  performed  in  the  five  subjects 
analyzed. 

D.  Data  Pre-Processing. 

Tirre  varying  spectrograms  of  either  full  and  low  resolution 
SL-transformed  EEG  data  by  estirrating  the  Power  Spectral 
Density  (PSD)  of  2-second  long  epochs,  each  starting  Is 
after  the  previous  one  were  computed.  The  Welch 
periodogram  algorithm  to  estimate  the  PSD  was  applied. 
Epochs  are  divided  into  segments  of  1  s,  with  a  Flann  window 
of  the  sarre  length  applied  to  each  segrrent,  and  50  % 
overlapping  between  the  segments.  This  gives  a  frequency 
resolution  of  1  FIz.  Finally,  the  power  components  are 
referred  to  the  correspond! ng  values  of  the  esti mated  PSD  of 
the  baseline  and  transforrred  in  dBo  i.e,  we  take  the 
logarithm  of  the  division.  The  spectral  values  were 
considered  i  n  a  frequency  band  from  8  to  30  FIz,  si  rce  those 
band  was  recognized  to  be  useful  for  the  recognition  of 
mental  pattern  i  n  pre/i  ous  papers  [  15, 16] . 

E.  Signal  Space  Projection 

I  n  the  Signal  Space  Projection  rrethod  a  w-dirrensional  space 
is  (Mined  so  that  a  i  measura  vector  (t),  whose 
components  are  features  extracted  from  incoming  data,  is 
represented  in  that  space  by  a  point.  I n  the  present  case,  the 
measure  vector  (t)  is  the  SL-transforrred  spectral  EEG  data 
in  the  frequency  band  of  8-30  FIz  computing  during  the 
mental  tasks  analyzed.  Given  p  vectors  of  w-dirrensional 


I  patterna  (  i,  2,0 ,  p),  the/?  components  of  the!  activation! 
vector: 

(t)=  (t)  (1) 

weight  the  preserce  of  each  pattern  in  (t).  is  the 
pseuctoi  nverse  of  the  proj  ecti  on  matri  x  whose  col  umns  are 
the  patterns  (  i,  2,0  ,  p).  The  pattern  descri  bi  ng  one  of  the  i**^ 
©cperimental  tasks  i  is  the  mean  of  the  selected  components 
of  the  PSD  computed  while  subject  was  irragining  or 
perforiri  ng  the  correspond!  ng  si  ngle  hand  rrDvement  (right  or 
left)  (training  procedure). 

F.  Fisher’s  linear  discriminant 

The  sarre  input  pattern  array  {t)  described  before  can  be 
cl  assif i  ed  with  a  general  I  i  near  di  scri  rri  rent  f  uncti  on  such  as 

y{t)=  '  it)  (2) 

where  is  the  array  of  unknown  weights  that  (Mines  the 
separation  between  classes  of  right  imagined  movement  (Rl) 
arid  left  irregined  moverrent  (LI)  in  the  input  space  It  is 
possible  to  (Mine  a  projection  that  maximizes  the  separation 
between  the  classes.  Fisheris  discriminant  [9, 10]  maximizes  a 
function  j(  )  that  represents  the  differences  between  the 
projected  class  rreans,  norrralized  by  a  measure  of  the 
within-class  scatter  along  the  w  direction.  By  Mining  r  and 
L  as  the  class  rreans  composed  by  the  sarre  spatial  patterns 
for  the  application  of  the  SSP  described  before,  and  the 
average  of  the  data  sets,  the  function  proposed  j(  )  is 

J(  )  =  ^ -  (3) 


where  i  s  the  between  cl  ass  covari  ance  rratri  x  gi  ven  by 

b  ~  L  fl  ■  L  fl 

and  Rw  the  withi  n-cl  ass  covari  ance  rratrix  given  by 
w=I(  it)-  ,)■{  it)-  J 

(5) 

+  Z(  it)-  l)  (  it)-  J 

leJi 

where  the  first  summation  runs  over  all  the  patterns 
belonging  to  the  class  describing  the  right  mental  imagery 
and  the  second  over  all  the  patterns  belonging  to  the  class 
related  to  the  Mt  rrental  imagery,  j  (  )  is  maximized  when 
the  weights  are  chosen  proportionally  to  the  -  r)- 

Oixe  obtained  the  weight  vector  w,  the  generic  input  pattern 
(0  wae  assigned  to  the  class  regarding  the  Mt  mental 
irragery  if  the  result  of  the  projection  '•  (?)-  is 

greater  than  zero,  and  to  the  class  related  to  the  right  mental 
irragery  otherwise 


G.  Classification  of  mental  patterns 

After  PSD  values  wae  computed  on  full  and  low  resolution 
SL-transforrred  EEG  data  such  values  are  then  fed  to  into 
the  linear  clasifiers  used  in  this  paper,  based  on  the  Signal 
Space  Projection  (SSP)  and  Fisher  discriminant  technique. 
These  dasifiers  wae  then  used  for  the  separation  of  rrental 
pattans  related  to  the  imagination  of  right  (Rl)  and  left  (LI) 
hand  moverrents. 

H.  Cross  validation. 

For  recognition  purposes,  we  applied  to  the  low  and  full 
resolution  SL  EEG  data  the  k-fold  cross-validation,  with  k  = 
8.  Flence,  we  divided  the  EEG  data  set  for  each  subject  i  nto  k 
subsets  of  equal  size.  The  SSP  and  Fisha  linear  discriminant 
projection  wae  raomputed  k  times,  each  time  leaving  out 
one  of  the  EEG  data  subsets  from  the  training,  and  using  the 
orritted  subset  to  corrpute  the  recognition  rate  Then,  the 
results  presented  hae  ae  an  avaage  of  the  recognition  rate 
obtained  for  each  one  of  the  k  subset  of  EEG  data  not  used 
for  the  SSP  and  Fisha  estirretion  of  the  das  means 
(training). 

I.  Statistical  analysis. 

A  two  way  Analysis  of  Vaiarce  (ANOVA)  was  performed 
on  the  avaage  values  of  the  recognition  scores  obtained  by 
the  cros-validation  technique  The  first  rrein  factor  was 
METFIODS  with  two  lei/els(SSP  and  FISFIER)  forthelinear 
classifiers  used  in  the  present  work,  while  the  saond  main 
factor  was  SPATIAL  FILTER,  with  two  levels  (LOWRESL 
and  FIIGFIRESL)  in  which  the  diffaent  implementation  of 
the  corrputation  of  the  surface  Laplacian  ae  corrpared  (low 
and  high  resolution  SL).  No  sphaical  correction  has  been 
used  [17]  fortheANOVA  corrputation  si  nee  the  lei/els  of  the 
mai  n  factors  ae  I  ess  than  three 

III.  Results 

Table  1  reports  the  avaage  correlation  values  obtained 
between  the  unprocessed  EEG  wavd'orms,  the  low-resolution 
and  full -resolution  SL  transforrration  of  the  EEG  waveforms, 
in  eah  time  poirt  aquired  and  on  all  the  five  raorded 
subjects.  Corrdation  values  wae  corrputed  for  the  channels 
included  in  the  computation  of  the  low- resolution  SL,  narrely 
F3,  Fz,  F4,  C3,  Cz,  C4,  P3,  Pz,  P4.  The  first  row  of  the  Table 
1  shows  the  correlation  values  between  the  low-resolution  SL 
and  the  full -resolution  SL  on  each  channel  analyzed.  Avaage 
correlation  between  low  ard  full  resolution  SL-transfomred 
EEG  wavdorms  was  0.65,  while  was  0.36  between  the 

Table i 

FIRST  ROW:  CORRELATION  VALUES  BETWEEN  THE  LOW-RESOLUTION  SL  AND  THE 
FULL-RESOLUTION  SL  (SL9-SL).  SECOND  ROW:  CORRELATION  VALUES  BETWEEN  THE 
LOW-RESOLUTION  SL  AND  THE  RAW  POTENTIALS  (SL9-P).  THIRD  ROW:  CORRELATION 
VALUES  BETWEEN  THE  SL  ESTIMATION  OBTAINED  USING  ALL  26  AVAILABLE 
CHANNELS  AND  RAW  POTENTIALS  (SL-P) 
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unprocessed  EEG  and  the  low-resolution  SL  wavdorms,  and 
was  0.38  between  the  unprocesed  EEG  and  the  full 
resolution  SL  wavdorms.  Table  2  reports  the  recognition 
scores  (in  pacentages)  of  the  mental  imagination  of 
movements  in  the  five  subjects  analyzed  with  both  the  linear 
classifiers  used  (SSP  and  Fisha)  with  data  from  low  and  full- 
resolution  SL.  With  SSP  using  the  low  resolution  SL  EEG 
data  the  avaage  raognition  was  of  81.3%  while  using  the 
full  resolution  SL  improves  to  82.1%.  The  use  of  Fisha 
classifia  applied  to  the  low  resolution  SL  EEG  data  produces 
60%  or  recognition  scor^  while  when  the  full  resolution  SL 
data  was  used  this  percentage  ari  ves  to  70.4%.  The  ANOVA 
demonstrated  that  the  use  of  SSP  improves  significantiy  the 
recognition  score  with  respect  the  use  of  Fisha  discriminant 
method  (METFIODS  main  factor,  F  =  11.75,  p<0.0266). 
Instead,  the  use  of  the  low  resolution  SL  data  does  rot 
decrease  significantly  the  raognition  rate  of  the  mertal 
patterns  with  respect  to  the  use  of  full  resolution  SL 
(SPATIAL  FILTERS  main  factor,  F  =  3.79,  p  =  0.12). 
Furthamore,  no  intaaction  METFIODS  x  SPATIAL 
FILTERS  was  found  (F  =2.90,  p  =0.16). 

Table  ii 

RECOGNITION  SCORES  FOR  THE  DETECTION  OF  RIGHT  AND  LEFT  IMAGINED 
MOVEMENTS  IN  FIVE  SUBJECTS.  PERCENTAGES  ARE  OBTAINED  WITH  THE  USE  OF  THE 
SSP  LINEAR  CLASSIFIERS  WITH  THE  LOW  RESOLUTION  SL-TRANSFORMED  EEG  DATA 
(SSP  LOWRES  SL)  AND  WITH  THE  FULL  RESOLUTION  SL-TRANSFORMED  EEG  DATA 
(SSP  FULLRES  SL),  AS  WELL  AS  WITH  THE  FISHER  LINEAR  CLASSIFIER  ON  LOW  AND 
FULL  RESOLUTION  SL-TRANSFORMED  EEG  DATA  (FISHER  LOWRES  SL  AND  FISHER 
FULLRES  SL,  RESPECTIVELY ) 
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83% 
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89% 

92% 
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85% 

Mean 

81% 

82% 

60% 

70% 

IV.  Discussion 

The  results  of  this  study  suggested  that  i  n  the  BCI  f rarre/vork 
it  is  usdul  to  corrpute  the  surfae  Laplaian  by  sphaical 
spline  also  if  a  limited  numba  of  scalp  electrodes  are 
avai lable  or  used  for  the  analysis  I  n  fact,  the  acuray  of  the 
computed  I  ow-resol  uti  on  SL  seems  to  be  not  too  far  from  that 
of  the  SL  obtai ned  by  26  seal p  electrodes  (avaage  correl ation 
coefficient  0.65  on  all  the  subjects  and  on  all  the  time  points 
analyzed).  It  is  rdei/ant  that  the  SL-transforrred  wavdorms 
with  both  rrodalities  (low  and  full  resolution)  showed  a  very 
low  correlation  with  the  unprocessed  raw  potentials  (about 
0.38  for  all  SL-transforrred  potentials).  More  important) y,  the 
avaage  mental  pattans  raognition  score  ova  five  subjects 
and  for  the  SSP  classifia  obtained  with  the  use  of  low- 
resolution  SL  ae  close  to  that  corrputed  with  the  SL 
computed  from  all  the  26  electrodes  used  (81.3%  and  82.1%, 
respotively).  Statistical  analysis  performed  with  the 
ANOVA  demonstrated  that  the  use  of  low  resolution  SL  data 
does  not  darease  significantly  the  performarce  of  the 
particula  dasifiers  used  (i.e  SSP  or  Fisha,  p  =0.16).  This 
result  is  promising  for  the  realization  of  BCI  devices  that  rely 


on  the  use  of  a  lirrited  number  of  electrodes,  also  at  the 
expenses  of  a  minor  recognition  rates  that  rraybe  can  be 
compensated  by  a  relative  larger  training  of  the  experimental 
subject.  Such  de/ices  can  compute  the  SL  transformed  EEG 
data  by  using  global  spline  SL  techniques  with  a  modest  loss 
of  accuracy  in  recognition  rates  of  mental  patterns  with 
respect  to  the  case  of  in  which  a  large  array  of  electrodes  is 
required.  Statistical  analysis  also  suggests  the  superiority  of 
the  Signal  Space  Projection  as  a  method  for  the  detection  of 
mental  patterns  with  respect  the  other  linear  discriminant 
technique,  namely  the  Fisher  linear  discriminant.  Compared 
to  neural  networks  [10,12],  I  inear  classifiers  are  easier  to  train 
since  they  do  not  require  non-linear  minimization.  With 
respect  the  recognition  scores  obtained  here,  other  Authors 
have  been  able  to  perform  successful  recognition  scores  of 
pattans  associated  with  the  preparation  of  performed 
movemants  with  linear  classification  technique  based  on  the 
Common  Spatial  Patterns  as  high  as  90%  (CSP)  [15]  as  well 
as  non  linear  classifiers  as  high  as  84%  [11-13]  in  the  BCI 
area 

I  n  surntTBry,  results  of  the  present  work  suggest  that  a 
BCI  device  based  on  linear  classifiers  and  Laplacian- 
transformed  EEG  signals  computed  froma  limited  number  of 
scalp  electrodes  (nine)  can  be  able  to  detect  mertal  activity 
with  a  reasonable  level  of  percentage  score  This  open  the 
avenue  for  more  practical  BCI  de/ices  that  does  not  requires 
the  use  of  a  I  arge  set  of  electrodes. 
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